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Abstract 

The dependency of the mehing point of supported metal nanoclusters as function of 
clusters height is theoretically investigated in the framework of the uniform approach. The 
H I vacancy mechanism describing the melting point shift in nanoclusters with decrease of their 

c/3 i size is proposed. It is shown that the essential role in clusters melting point shift is played 

by van der Waals forces of cluster-substrate interaction. It is shown, that the account of 
layer-by-layer fusion of a cluster allows to satisfactorily describe the melting of nanoclusters 
of various metals, deposited onto a different substrates. The proposed model satisfactorily 
accounts for the experimental data. 
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O , One of the interesting properties of nanometer-sized systems is a shift of the melting point 

and the lattice parameter with decrease of cluster size (size effect). This effect was experimentally 
_ observed for Au, Ag and Cu (THH] clusters deposited onto various substrates (W, C, Ni). The 

^ ! 

t~^ ! existing theoretical descriptions of the observed size effect are based on various approaches [H 

[■^ , 

l/^ . [2|[7]-|10]. However, their common and essential drawback is that the influence of the substrate 

is not taken into account. At the same time it is experimentally proved [6] that substrate plays 
O I a crucial role and the use of different substrates (W, C) results in the significant difference in 

melting points of clusters (up to 20% for 20 A clusters on carbon and tungsten substrates [6]). 

Physically it can be related to the additional energy of the supported cluster due to the contact 
/\^ • potential difference [llj, the lattice mismatch at the cluster-substrate interface [T2|, and the 
C^ ' cluster-substrate van der Waals interaction [T31 that can lead to the shifts of cluster's lattice 



parameter and melting temperature. Therefore, the understanding of the physical reason for the 
size shift of cluster's melting point and lattice parameter needs the determination of the influence 
of the substrate on cluster's properties. 

In this paper the theoretical description of the melting point and the lattice parameter shifts 
for the supported nanoclusters of various metals carried out in the framework of the uniform 
approach using the vacancy mechanism taking into account the cluster-substrate interaction is 
proposed. In the framework of this mechanism the change of cluster properties with decrease of 
its size is described as a result of anharmonic oscillations of cluster atoms, leading to the decrease 
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of the energy of vacancy formation and, hence, to the additional creation of vacancies in cluster 
with decrease of its size. 

Recently in this paper [T3] the theoretical description of the melting point for the supported 
nanoclusters of various metals carried out in the framework of the uniform approach using the 
vacancy mechanism taking into account the cluster-substrate interaction is proposed. In the 
framework of this mechanism the change of cluster properties with decrease of its size is described 
as a result of anharmonic oscillations of cluster atoms, leading to the decrease of the energy of 
vacancy formation and, hence, to the additional creation of vacancies in cluster with decrease of 
its size. 

In the paper [13] the qualitatively agreement with experimental data of melting of Au clusters, 
deposited onto various substrates (W, C) is obtained. The difference may be related to the 
cluster-substrate lattice mismatch, to the contact potential difference at the interface, as well as 
to the fact that melting of the cluster doesn't occur in the whole volume of cluster but begins 
with the surface layers while the inner layers remain solid (that means the existing both the 
liquid surface layer and solid core, "layer-by-layer melting" below). The letter effect is observed 
experimentally [6],|T5[|1^. Thus, the quantitative description of the melting temperature shift 
needs the determination of the influence of these effects on cluster's properties. 

In this paper the influence of the contact potential difference at the interface, the finite lateral 
cluster size and the layer-by-layer melting of the clusters had been taken into account so that we 
succeeded in quantitative description of the melting temperature shift of different metal clusters 
deposited onto different substrates (Au, Ag clusters deposited onto surfaces of W, HOPG, Ni). 

This effect can be described considering the van der Waals interaction of a supported metal 
cluster with bulk metal or semimetal substrate |13j . Van der Waals forces are of the electro- 
magnetic character and arise due to the mutual polarization of neutral atoms of the interacting 
objects. The presence of these forces can lead to the change of the chemical potential of a 
cluster. In this case, if the change is positive, it can be interpreted as an additional thermody- 
namic potential per atom of a cluster, arising due to the interaction with substrate. It results 
in the decrease of the binding energy of atoms in a lattice and, therefore, in the increase of the 
probability of vacancies formation in a cluster. Besides van der Waals forces, the formation of 
vacancies in a cluster occurs due to the presence of cluster surface, that results in the dependence 
of the energy of vacancy formation on its distance from cluster surface, as well as on cluster 
shape [17J. Under assumption that the clusters under investigation jIHS] satisfy the requirement 
(where / = 0.7 — 10.0 nm and h = 0.3 — 5.0 nni are cluster's lateral size and height, consequently), 
the cluster surface can be considered as a fiat and the influence of the boundary effects related 
to cluster's shape can be neglected. 

The equilibrium concentration of vacancies in a cluster can be found from the requirement 
of the minimum of the free energy of cluster's vacancies subsystem taking into account the 
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interaction with substrate: 

71=71 

where n is the average concentration of vacancies in a cluster. The interaction of vacancies at 
the distance more than the character size of the vacancy (the possibihty of divacancy formation) is 
attractive, while at the distance less than the characteristic size of the vacancy it has a character 
of a solid core (it is impossible to locate two vacancies at the same site of the lattice). Assuming 
the equilibrium concentration of vacancies much less than the atomic density in a cluster, the 
vacancies subsystem of a cluster can be considered as a van der Waals' gas [T8]. In this case the 
free energy of the vacancies subsystem after decomposition for virial coefficients is given by |19] : 

,.,^,!flin, B(r)4/(l-exp(-^)).V 

where N is the number of vacancies in a cluster, V is the cluster volume, U12 is the interaction 
potential of vacancies, and Fid is the free energy of non- inter acting vacancies system. Considering 
U12 as a rectangular potential well of depth Ed (the energy of divacancy formation), the virial 
coefficient B{T) is expressed as: 

S(T) = 6 + 76(l-exp(E,/T)) (3) 

where h = —p^, tq — is the characteristic size of a vacancy. The free energy Fid of non-interacting 
vacancies system is given by |T7] : 

Fid = / EyU dV + T / [n In u/uq — n] dV 

J J (4) 

E, = Ef - 6fi^''^ + 6fi, 

where Eb, Es are the energies of vacancy formation in the volume and at the surface of a 
cluster, consequently, tiq is the atomic density for a cluster, a is the cluster lattice constant, 5^'"'^^ 
is the additive chemical potential of a cluster arising due to its interaction with substrate |13| . 5/ic 
— the shift of the cluster chemical potential due to the presence of the cluster-substrate interface. 

The energy of the vacancy formation E^ is assumed depending on its distance from cluster 
surface [T7] : 

Ef = Eb- {Eb - Es) exp (-^) (5) 

where Eb-, Es are the energies of vacancy formation in the volume and at the surface of a cluster, 
consequently, h is the cluster height, a is the cluster lattice constant, x is the distance from the 
cluster surface. Substituting the expression (jlj) and into ([2]), the free energy is obtained: 



F = E^ndV + T [nlnn/no -n] dV - nTln (1 - nb) - 7n%T (1 - exp (Ed/T)) (6) 

The equation for the equihbrium concentration of vacancies rl in a cluster follows from ([6]) 
and ([T]) as: 

^B _ g^vdv + 5^^ + 2- In iL _ T In (1 - nb) + nT—^ - UnbT [1 - exp ( ^] ) = (7) 

no 1-nb \ \T J J 

where Six"'^^ can be calculated if the complex dielectric constants of substrate and cluster are 
known fT3]. The effect of the cluster surface taken into account in the expression for the energy 
of vacancy formation (jl]), the value 5^'"'^'" can be calculated for a film of thickness h deposited 
onto substrate [13]: 

A„-^J^. rr- [ {ef{iO-l){ef{ii)-es{iO) . 

^ zSn^h^' "^ J {efm + l){efm+esm) ^ ^^ 



where i^ = u is the imaginary part of the electromagnetic field, Eg, Sf are the dielectric 
constants of substrate and cluster materials, consequently, and z is the number of the nearest 
neighbors in a cluster. Under assumption of the validity of the expression Ej (z^) = 1 + AnkY, 
j = s, f for the dielectric constants [I7l|20], k = 8.98 • 10^ A^ ■ m? /C'^ the value dfi"'^'" is given by: 



.vdv _ «^ ^(^/ - ^s) (Tf . f(^f + (^s 



^ vctv ^ _ -"v-/ -^^ . r/ . In Il^^ll (9) 
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where o"/, (Jg are the conductivity of cluster and substrate, consequently. 

The shift of the cluster chemical potential due to the presence of the cluster-substrate interface 

is given by |11] : 

5iic = ^exp {—K{h — x)) 

_ K2 (/i2 - /il) 

K + /.2) (10) 

n?2 

K?n = 47re^— 4- 

^'2 kT 

where fii, ^2 are the chemical potentials of graphite (1) and Au (2), n\ I is the electron density. 

The conductivity o"/ of a cluster can be calculated by taking into account the contributions to 
the intrinsic resistivity of the cluster material po and the additional resistivity p^ac arising due to 
the additional scattering of conduction electrons from vacancies. The additional resistivity p^^c 
depends on the concentration of vacancies defined by the additional chemical potential arising 
due to the presence of a substrate. Thus, the expression for the cluster conductivity takes the 
form: 
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^/ = (PO + Pvac) ; Pvac = pUc " ^/^O 

where pmc is the additional resistivity per a vacancy that does not depend on n. The ex- 
pressions (E]), ([9]), (ITTj) allow to calculate the equilibrium concentration of vacancies in a cluster 
taking into account the cluster-substrate interaction. The theoretical dependence of the relative 
concentration of vacancies for Au cluster on graphite surface as function of cluster height at room 
temperature (T=300°-ft') calculated from the equations (jTj), ([9]) and (fTTj) is plotted in Fig. [1]. 
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Figure 1: The equilibrium concentration of vacancies n in Au cluster on graphite at T = 300° X relative 
to the atomic density uq calculated as function of cluster height h. 




It is seen that the concentration of vacancies sharply increases for cluster height less than 22A, 
that can result in the modification of the modulus of elasticity of cluster and, as a consequence, 
to the shifts of the lattice parameter and the melting point [18j. 

The melting point can be calculated taking into account the dependence of the modulus 
of elasticity on the concentration of vacancies in a cluster [T8]. Using the vacancy model of 
melting [18] and considering vacancies as dilatation centers, it is possible to calculate the shear 
modulus as a function of the concentration of vacancies [18j: 
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K - A.^ {K/^Vf - n{l-n/n,f 

15 r[l + 7n/no(l-n/no)^ 



;i-exp(E,/r))] 
Here K is the shear modulus of the defect free material and AV is the dilatation volume. 
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Taking into account the finite lateral cluster size, the dilatation volume depending on the 
particle size can be found. The stress tensor cr*^ in an elastic medium with point defects can be 
calculated in the model considering the defect as centre of a dilatation [21j : 



"•i = A' - -f, S„S(r) + 2f,^ - 6f,iiii AV (13) 



where /i is the modulus of dilatation. 

On the other hand, let's consider a spherical cluster with point defect as a hollow ball with 
exterior radius i?i, interior radius R2 with pressure outside pi, and pressure in a cavity P2 = 0. 
In this case strain of a hollow ball is given by |22] : 



where a = ^ 



W^2 



Ui= ( a + — j Xj, (14) 

and h = ^-nrzm has been found from boundary conditions [22j. In this 



case the stress tensor aik takes the form: 

a., = {{l{^- l^') ^ (^ + 3A^«) S^k + 2/i^ - 6/i^} b (15) 

Equations ( IT3l) and (ITSl) allows to calculate the dilatational volume AV as a function of the 
cluster radius: 



A^ = ^=f^#^ (16) 



P2 UlRf 
4/i i?3 - R^ 

The dilatational volume AVq for a continuous medium with a radius of cavity i?2 can be found 
from boundary conditions (i?i -^ 00, AV -^ AVq): 

AK, = T^i?^ , (17) 

4/i 

Finally the dilatation volume AV^ takes the form: 

Substituting the expression fITS]) into 012p the shear modulus is obtained: 

A- = A- - ±.' (k^V„^-^) ' -^ ni^^nln^ ,jg, 

15 V Ri-RlJ T[l + 7n/no(l-n/naf(l-exp{EJT))] 

The dependence of the shear modulus K of Au clusters supported on carbon sas a function of 
cluster height L at temperature T = 300° is presented in Fig. [2l Defining the melting point of a 
cluster from the Born's criterion K{Tmeit) = jl8], the dependence of the melting temperature of 
a layer at a distance x from a cluster surface on cluster height can be obtained from the equation 



(1121) . The cluster is considered to be melted if it has no solid core. The dependencies of the 
melting temperature of the last layer (the nearest to the cluster - a substrate interface) for Au 
and Ag clusters deposited onto various substrates (HOPG, W, Ni) as function of cluster height, 
calculated from ( !T9|) are presented in Fig. |3l 
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Figure 2: The dependence of a shear modulus K /K^ where X is a shear modulus of defect-free substance 
on the cluster height H at temperature T = 300° ET 
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Figure 3: The dependence of the melting point Tmelt of different clusters on its height H . 
are the theoretical dependences calculated using equations p9|) for Au, Ag clusters on W, 
Points are the experimental data (Au/HOPG, Au/Ni [23], Au/W, Ag/W [6j) 



Solid curves 
HOPG, Ni. 



It is seen from Fig. [3] that an additional formation of vacancies in a cluster due to the inter- 
action with substrate, leads to reduction of melting temperature that quantitatively agrees with 
the experimental data. 

Let's consider the dependence of melting temperature of a cluster on distance to a cluster 
surface. The theoretical dependence of melting temperature of layers of Au cluster oi H = 60Ain 
height, deposited onto W surface, as a function of a distance to the cluster surface is presented 
in fig. m 
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Figure 4: The dependence of melting temperature of the layer which as a function of distance x to the 
surface of Au cluster of H = 60 A in height, deposited onto W surface 

It is seen from fig. HI that melting temperature of surface layers is considerably lower (^1.5 
times), than melting temperature of the layers neighboring to a cluster - substrate interface. 
Surface of a cluster of 60A in height is considered to be melted already at T = 800° K whereas 
the inner layers are melted at considerably higher temperatures, so that the melting of a whole 
cluster with height of 60A occurs at temperature T ^ 1200°i^. This behavior is observed exper- 
imentally for various materials of clusters (see for example [61 IT51IT6] ). 



Let's estimate contributions of the considered effects to the melting temperature of clusters: 
the additional energy of a cluster per one atom due to its interaction with substrate can be 
estimated in the frames of several mechanisms taking into account the contact potential difference, 
the lattice mismatch at the interface, the van der Waals interaction with the substrate, as well 
as the layer-by-layer melting. 

The estimations will be carried out for the cluster of 30A in height. According to the equations 
([9]) and (ITTl) , the additional energy of a cluster per one atom due to the van der Waals interaction 

with substrate is (^/id^Lr L=3oA ~ ^■'^^^■ 

The additional energy of a cluster obtained due to the contact potential difference can be 
estimated assuming that the influence of a substrate is considerable only within the flrst monolayer 
of atoms at the interface: 5iXc ~ ^ / Sfj.dV ~ ^^Jj.'^h ~ h^ ' "^^ere V ^ h are the volume and the 
height of a cluster. In this case the value of the additional energy due to the contact potential 
difference equals to ^^Kiuster ~ — O.OSeV. 

/i=30A 

The additional energy of a cluster per one atom due to the lattice mismatch at the interface 
can be estimated considering that the characteristic temperature of the transition to the dispro- 
portionate phase is about ^ bOO°K ^24j. Thus, the additional energy per one atom due to the 



this mechanism is 5^^°j^^lter ~ 0.05e\^. 

In this paper the energy of vacancy formation had been taken into account to be different for 
a surface and bulk, in contrast to the previous paper jl3] where the energy of vacancy formation 
was assumed to be mean quantity of a cluster height. Thus, according to the equation ([5]) 
the additional energy of a cluster of 30Ain height due to the different formation energy of the 
vacancies for the surface and bulk is 5E\^^^^y^ ^ {].leV . 
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Figure 5: The dependences of the melting temperature of Au cluster deposited onto HOPG surface. A 
solid curve — with the account of the layer-by-layer melting, a dotted line is the dependence calculated 
without the account of the layer-by-layer melting, a dashed curve is the dependence calculated without 
the influence of a contact potential difference. Points — are the experimental data |23] 

Fig. |5] shows the dependences of melting temperature of Au clusters deposited onto HOPG 
surface on the cluster height taking into account the layer-by-layer melting (a continuous line), 
without taking into account layer-by-layer melting (averaging energy of vacancy formation on 
cluster height) (a frequent dotted line) and without contact potential difference (a long dotted 
line). It is seen that taking into account the layer-by-layer melting results in significant changes 
in the curve behavior that quantitatively agrees with the experimental data. Taking into account 
the contact potential difference of a cluster - substrate changes a trend of curve insignificantly. 
Taking into account the finite lateral size of a cluster had not given appreciable changes in melting 
temperature of a cluster. 
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Thus, it is shown, the major effect on the shift of the melting point and the lattice parameter 
of supported nanoclusters under the high vacuum conditions is produced by the van der Waals 
cluster-substrate interaction. Taking into account the contact potential difference on a cluster- 
substrate interface, finite lateral cluster size, and layer-by-layer melting of nanoparticles has 
been considered allowed one to satisfactorily describe the dependences of melting temperature 
on cluster height for Au and Ag clusters, deposited onto surfaces of W, HOPG, Ni. 

It is shown, that the account of layer-by-layer fusion of a cluster (i.e. observationally found 
phenomena when the melting of a cluster doesn't occur simultaneously in the whole volume 
of a cluster, but the presence both a liquid surface layer and a solid core is possible) gives the 
considerable contribution to melting temperature of the clusters deposited onto various substrates 
that has allowed to one to satisfactorily describe the melting of nanoclusters of various metals, 
deposited onto a different substrates. 

This work was supported by the Federal target program "Scientific and scientific-educational 
personnel of innovative Russia". 
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